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Kenneth J. Baumeister 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

A finite-element Galerkin formulation has been developed to study trans- 
verse magnetic (TM) wave propagation in two-dimensional S-curved ducts with 
both perfectly conducting and absorbing walls. The reflection and transmission 
at the entrances and the exits of the curved ducts are determined by coupling 
the finite-element solutions in the curved ducts to the eigenfunctions of an 
infinite, uniform, perfectly conducting duct. Example solutions are presented 
for a double mitred and S-ducts of various lengths. The length of the S-duct 
is found to significantly effect the reflective characteristics of the duct. 
Also, the effect of curvature on an absorbing duct wall is illustrated. 


INTRODUCTION 

Electromagnetic propagation in curved ducts (wave guides) plays an impor- 
tant role in many practical physical systems. In microwave power generation 
systems for example (ref. 1, pp. 315-317), bends or corners are required to 
alter the direction of the wave. Since curves or corners represent discontinu- 
ities, reflection from the bends can be significant. To better understand the 
electromagnetic transmission properties of a bend, a finite-element Galerkin 
formulation has been developed to study wave propagation in curved S-shaped 
ducts. Both perfectly conducting and absorbing walls will be considered. 
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mode amplitude of positive going entrance waves, equation (15) 

mode amplitude of reflected negative going entrance waves, 
equation (15) 

mode amplitude of positive going exit waves, equation (18) 
characteristic duct height 

dimensionless entrance height b'/b' 

dimensionless exit height b^/b ' 
speed of light in vacuum 

dimensionless harmonic electric field vector, 
unit vectors in coordinate directions 
dimensionless frequency, equation (7) 
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dimensionless harmonic magnetic intensity vector, H'(x,y,z)/iV 

x component of magnetic intensity H 

magnetic intensity at node i 

normalizing magnitude of magnetic intensity 

finite-element approximation to H 

wave number, w/c 

axial modal wave number, equations (16) and (17) 

dimensionless length, L'/b' 

number of elements 

mode number, equation (20) 

number of nodes 

number of modes 

local interpolation shape functions 

mode number, equation (15) 

outward normal unit vector 

length of line segment on boundary 

c ' 

dimensionless time, t' 

weight, equation (25) 

dimensionless transverse distance, x'/b' 

dimensionless transverse distance, y'/b' 

dimensionless axial distance, z'/b' 

dielectric constant, e'/e^ 

dielectric constant in entrance duct, e'/e' 

a 0 

permi tti vi ty 

permittivity in vacuum 

complex permittivity, equation (5) 
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P r relative permeability, p'/p^ 

p 1 dimensional permeability 

p^ permeability in vacuum 

a dimensionless conductance, a'b'/c'e' 

o o 

co' angular velocity 

a) dimensionless angular velocity, w'b'/c^ 

Subscripts 

a entrance region 

b exit region 

x,y,z scalar vector components 

Superscripts 

(1) region 1 

(2) region 2 

' dimensional quantity 

T transpose 

vector quantity 


GEOMETRICAL MODEL 

In the finite-element modeling of the curved ducts to be presented, an 
S-shaped profile has been closed to approximate the two-dimensional cross- 
sectional profile that might be found in a typical bend, as shown in figure 1. 
The S-shaped profile can be prescribed by a simple third-degree polynomial of 
the form 


y ■ 3 U 


il 


(l) 


where the dimensionless duct coordinates are defined as 


y = 


. Il 2 = 


b a 


L - 


b a 


< 2 > 


and b' is the height of the straight entrance duct leading into the curved 


duct. The S-curve defined by equation (1) has zero slope at z/L of 0 and 1; 
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providing a smooth transition from a straight entrance to the curved test sec- 
tion. In the foregoing equations, the prime, ', is used to denote a dimensional 
quantity and the unprimed defines a dimensionless quantity. This convention 
will be used throughout this paper. These and all other symbols used in the 
report are defined in the nomenclature. 

This paper will focus on the interaction of a propagating duct mode travel- 
ing down the uniform entrance duct with the curved wall as shown in figure 2. 

The reflection and transmission at the entrance and exit of the curved duct are 
determined by coupling the finite-element solutions in the curved duct to the 
eigenfunctions of the infinite, uniform, perfectly conducting entrance and exit 
ducts. This permits a multimodal representation accounting for reflection and 
mode conversion by the nonuniformity (ref. 2). 


GOVERNING EQUATIONS 

The governing differential equations are the standard Maxwell's equations 
along with the boundary conditions associated with perfectly conducting and 
lossy walls. Maxwell's equations will be combined to form a single wave equa- 
tion for transverse magnetic wave propagation in a two-dimensional duct. 


Maxwel 1 ‘ s Equations 

Maxwell's curl equations for time harmonic (e + 3 ut ) variations in the elec- 
tromagnetic fields are 


VxE = -jcjprH 
VxH = jcoeE 

where the total permittivity included conduction 

e = e f - ja/u 


a = 


°' b i 


and the dimensionless frequency is defined as 
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f ' b ' 
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(5) 

( 6 ) 

(7) 


Heterogeneous Variable Property Wave Equation 

The number of dependent variables can be reduced by combining Maxwell's 
equations (3) and (4) into a single wave equation. It is desirable to develop 
a wave equation that could be used for varying media properties so that no 
special treatment of the interface between materials is required. That is. 
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the same equations apply in the duct and in the absorber region and only the 
material properties are changed. 

Rewriting equation (4), 


VxH - 
jwe 


(8) 


and taking the curl 


Vx 


VxH 

jcoe 


VxE 


(9) 


The constants ju are independent of space and can be pulled out the curl 
operator in equation (9); however, e must remain inside since e is now 
assumed a function of the spatial dimensions. 


J_ 



VxE 


Substituting equation (10) into equation (3) yields our heterogeneous 
governing wave equation. 


(10) 
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( 11 ) 


TM Variable Property Two-Dimensional Wave Equation 

Transverse magnetic (TM) waves will now be assumed to represent the input 
electromagnetic modes propagating down the entrance duct towards the curved 
S-duct. For this two-dimensional geometry, the assumption is made that only 
one component of the H vector will exist in the problem domain, that is, 


H(y,z) = H x (y,z)e x (12) 

where e x is a unit vector in the x direction (into the paper as_shown in 
fig. 2). The magnitude of the single x component of the vector H depends 
only on the two spatial dimensions y and z. Equation (12) is a valid 
solution in both the straight and curved duct since it identically satisfies 
Maxwell's divergence equations for magnetic fields. 

Substituting equation (12) into equation (11) and performing the required 
manipulation yields 
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(13) 
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In vector form, 


V • - VH + co 2 u H = 0 (14) 

C X I A 

Equation (14) represents the governing wave equation to be solved by 
finite-element theory. 


UNIFORM DUCT ANALYTICAL SOLUTIONS 

The analytical solutions of equation (14) for wave propagation in the 
uniform perfectly conducting duct having an anechoic entrance and exit will be 
employed to give the termination boundary condition for the finite-element 
region. The analytical solution of equation (14) for TM waves traveling 
between perfectly conducting parallel plates is given as (ref. 3, p. 458): 
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/T cos ' "- yj e J, 'zn^ ’ ^ A p cos - -- - -- yj e^zn" (15) 

n=l 


For the e^ wt time dependence used here, the A*e ^ k zn z term represents a 

— i k z 

wave propagating in the positive z direction while the A n e J zn term 
represents a wave moving in the negative z direction. 

The axial wave number k zn in equation (15) is 
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A similar solution exists at the exit, except only positive going waves 
are considered. 


xb 


N. 


n=l 


8 + cos (I 4-U« y e -j k zn z 
n b, J 


BOUNDARY CONDITIONS 


(18) 


A variety of boundary conditions will be used in the finite-element solu- 
tion of equation (14) for the model problem which is displayed in schematic 
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form in figure 2. Each of the required conditions will now be briefly 
discussed. 


Input Condition 

The analysis assumes a given number N m of propagating A* modes 

(eq. 15). These modes effectively set the level of the magnetic field in the 
finite-element region and can be viewed as the equivalent Dirichlet boundary 
conditions required for elliptic boundary value problem as defined by 
equation (14). 

The modal expression represented by equation (15) has been truncated to a 

total of N modes of the infinite number possible. Thus, a total of N 
m _ m 

unknown modal amplitudes Aj, A^, . . . A^ have been introduced. N m con- 

m " 

straint equations will be required to determine each of these unknown reflec- 
tion coefficients. The equations used to define these coefficients will now be 
i ntroduced . 


Continuity at Inlet and Exit 

The tangential component of an H field is continuous across an interface 
between two physically real media which are not perfect conductors (ref. 4, 
eq. (1.61)). Thus, the boundary condition becomes 

H xa = H x (z = ° ; 0 < y < b & ) (19) 

where H is the modal representation of magnetic field in the analytical 

X cl ^ 

inlet region given by equation (15) and represent the finite-element 
approximation for H at the interface. The hat over H implies an approxi- 

X X 

mate finite element numerical solution to be discussed in detail in a following 
section. 

At the inlet to the curved section, shown in figure 2, the H xa in the 
analytical region given by equation (15) must match the magnetic field defined 
by the finite-element nodal points along the boundary interface. Many possible 
matching methods can by employed for this boundary condition, such as point 
collocation, least squares, or weighted residuals. A weighted residual approach 
was used herein with the weighting function equal to the eigenfunctions for the 
uniform infinitely long duct with perfectly conducting walls; 



(N equations, m = 1,2,3, N m ) 
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Equation (20) represents N m separate equations; one for each coefficient 
defined in equation (15). The symbol m has been introduced as a dummy vari- 
able to make it distinct from the multiple n mode numbers that make up the 
H a analytical function. A similar equation exists at the exit. 

In addition to the tangential component of the magnetic field, the tangen- 
tial component of the electric field must also be continuous across the inter- 
face (ref. 3, eq. (7-52 ( a) ) ) . Using equation (4) to express the tangential 
electric field in terms of the magnetic field (ref. 5, eq. (7-4)) yields 



e 3z e 

& 



( 21 ) 


The weak form of the finite-element solution will be employed in the solu- 
tion of this problem. In this form, a contour integral term will be developed 
which will contain a natural boundary condition of the form VH • n where n 
represents the unit outward normal. Equation (21) can be expressed in this 
gradient form at the entrance as 

8H 

VH • n = — (inlet) (22) 

X e, dZ 


and at the exits 


- +e 8H xb 

VH • n = ^ (exit) 

X dZ 


(23) 


The sign change in equations (22) and (23) comes directly from the directional 
change of the unit outward normal n. 


Perfectly Conducting Wall Conditions 

At a perfectly conducting wall, the tangential component of the electric 
field vector is zero (ref. 3, eq. (7-52(a)) or ref. 4, eq. (1.69)). Again, 
using equation (4) to relate the electric field to the magnetic field (ref. 5, 
eq. (7-4)), the component of the gradient of the magnetic field normal to a 
perfectly conducting wall becomes 


VH x • n = 0 


(24) 


FINITE-ELEMENT THEORY 

The finite-element formulation of the electromagnetic wave equation is now 
generated by using the method of weighted residuals to obtain an integral form 
of the variable property wave equation over the whole (global) domain D shown 
in figure 3. 
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System Di scretization 


The continuous domain D is first divided into a number of discrete areas 
staked out by the nodal points as shown in figure 2. 


Global Weighted Residual Approach 

In the^cl assi cal weighted residual manner, the magnetic field intensity 
component H x (y,z) is curve fitted by expanding in terms of all the unknown 

nodal values* H .(y. ,z.) and a known series of basis (weight) functions, such 
that xi t i 

N 

H x (y,z> = J W.(y,z)H x . = CW]{H x } (25) 

i-1 


where the basis or weight functions W.(y,z) characterizes the spatial depend- 
ence of H (y,z) in terms of H . which represents the unknown value of the 

magnetic field intensity component at the i 1 nodal point in the global 
region. As before in equation (19), the hat over the H x (y,z) indicates that 

it is the approximate numerical solution to H x (y,z). The precise form of the 
known weight which will be employed in this analysis will be presented later 
in the section entitled Galerkin Approximation. 

Upon substituting the assumed magnetic field series, equation (25), into 
the governing wave equation at each nodal point, a distribution of errors 
results throughout the finite-element region due to the approximate nature of 
the assumed series. In accordance with the method of weighted residuals, the 
integral of the assumed basis function W^ and the error at each nodal point 
are forced to be zero (orthogonal) by letting 

W. ( v * l v ^x + cj2 Vi r ^x ) dydz = 0 (26) 

JJ D V ) 

(i = 1 ,2, . . . N equations) 

Thus, there are N separate equations (written in compact form); one equation 
for each of the N nodal H x i unknowns. 

By making use of the gradient vector identity of a scalar and a vector 
and Green's theorem in a plane (ref 
be converted to 


1 vw. • VH x - W. u 2 p r H 

J D 

(i = 1,2, . 


7, p. 79, eq. (4.7(b))) equation (26) can 


dydz - <{> W. - VH • n ds = 0 (27) 

llC X i 


. . N equations) 
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In effect, the second order differential equation has been reduced to a 
first order equation allowing the use of the weak formulation of the finite- 
element theory. 


Finite Element Approximation 


Equation (27) is valid over the entire domain D shown in figure 3 or any 
subdomain A e , as represented by the area of a small triangular element embed- 
ded in the region as depicted in figure 3. To begin the finite-element aspect 
of the weighted residual method, the domain D is assumed to be divided into 
M elements defined by N nodes, see figure 2. In this case, equation (27) 
can be written as 


e=l 


'Ae 


1 ~(e) 

— VW. • VH ve ' 

£ 1 X 

x e 


II 2 S (e > 
M 1“ ' J re H x 


dydz - <b 


W. — VH (e) 
l e x 
v e 


• n ds = 0 


( i = 1 ,2 . . . N equations) (28) 

Where the properties are now given a subscript e to indicate that they may 
vary from element to element. Although the area integration in equation (28) 
applies to an individual element, the line integral is still defined over the 
global surface area, and as such, can be treated independently from the local 

~(e) 

area integrations. Also, H has been replaced by H which expresses the 

X X 

magnetic field in terms of the three local element nodes associated with each 
triangular element as shown in figure 3. 


Galerkin Approximation 

In the method of weighted residuals, the weight Wi is assumed to be a 
known function which has the property of being unity at node i and identical 
to zero at all other nodes. The Galerkin approximation to the more general 
weighted residual approach assumes that W^ can be related to local shape 

( p ) 

function N. (x,y) which represents the variation of the field variable H 

1 X 

and its derivatives inside the element. 


The form of the local shape matrix [N] (e) depends on the type of element 
used. For the linear triangular element employed herein, the known value of 

N (e ^ can be found in texts on finite elements (ref. 9, p. 6). Thus, the mag- 
netic field H x (y,z) inside a particular element can be expressed in terms of 

the local shape function N.(y,z) as follows: 


H (y,z) 
x J ’ 


H^ e) (y,z) 


,,(e), , u (e) ,\u (e) »i ( e ) ^ _,,. u (e) 

N ] (y,z)H xl + N 2 (y,z)H x2 + N 3 (y,z)H x3 


N$ e) H*f > = [N] (e) {H } (e) 

J x J X 


3=1 


( 29 ) 
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where {H x }^ e ^ is the vector of nodal values of H x for a general element e 
with subscripts 1, 2, and 3 representing the nodal positions as shown in fig- 
ure 3. The subscripts 1, 2, and 3 take on the actual nodal number of figure 2 
when applied to a specific element in the finite-element region. 

(e) 

The weight W.(y,z) is now approximated by multiple values of N. 

W.(y.z) = N^ e) (y,z) (30) 

in all elements containing the i^h node. For all elements which do not con- 
tain the node i, the weight W^ is assumed zero not only at all other nodes 
(as required by the general definition of Wi ) but also at all values of y 
and z in the elements which do not contain the node i. Thus, W^ can be 
visualized as a multisided pyramid on a flat surface with its apex over node i 
and corners at adjacent nodes. 

(e) 

Recognizing that N. is zero for all elements not having the unknown 
H . associated with a particular element, the finite element equation (28) can 
now be written in compact form as 


M 


e=l 



f— VN^ e) • VH <e> 
ve 1 


(e) 


2 

w (i 



dydz 


- d> ^lj e) ] — VH^ e) • nj ds = 0 (31) 


(i = 1 ,2, . . . N equations) 


Global Matrix 

The finite-element aspects of converting equation (31), the entrance con- 
dition (20), and the similar exit condition into a set of global difference 
equations can be found in text books as well as reference 10 and for concise- 
ness will not be presented herein. However, equation (70) of reference 10 is 
identical to equation (31). In reference 10, additional details can be found 
for the manipulation of equation (31) into the final finite-element solution. 


RESULTS AND COMPARISONS 

For theory and code validation, the finite-element solution is first 
applied to a straight walled case where a simple exact solution exists. Next, 
the solution is applied to a mitred duct approximation to an S-duct for which 
certain results can be expected from previous numerical solutions. As a third 
example, the finite-element solutions will be presented for a set of S-ducts of 
various lengths with perfectly conducting boundary conditions. Finally, the 
solution for a curved duct with an absorbing wall will be presented. With the 
exception of the first example, all ducts are air filled (e = e r = 1 , y r = 1). 
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Example 1: Reflection and Transmission with Normal Incidence 

Consider the case of an incident plane magnetic wave H x of dimensionless 
frequency 2ir encountering a step change in the dielectric constant from 1 to 
4 at the dimensionless axial position of z equal to 0.25 inside the finite- 
element grid. As shown in figure 4, the finite element and exact analytical 
(ref. 5, chapter 5) theories are in good agreement for the absolute magnitude 
of the magnetic intensity H x . Additional validation examples for propagation 
with straight walls can be found in reference 10. 


Example 2: Double Mitre Approximation to S-Duct 

As shown in figure 5, a double 90° mitred duct can be represented as a 
rough approximation to an S-duct. Plane wave acoustic propagation in a single 
90° mitred duct is well understood from the numerical studies of Shepherd and 
Cabelli (ref. 11). The energy reflection coefficient of a single 90° mitred 
bend for plane wave incidence is shown in figure 6 from reference 11. As the 
frequency of the incoming plane acoustic wave approaches the first mode cut-on 
frequency, the transmitted mode shape in the exit portion of the mitred duct 
approaches the shape of the first higher order mode. Consequently, as shown 
in figure 6, nearly all the incoming energy is reflected back down the duct, 
since the transmitted wave is cut-off and cannot carry energy. For very low 
frequency nearly all its energy is transmitted down the duct while at higher 
frequencies complete reflection of energy is expected, as depicted by the ray 
approximation shown in figure 6. 

As seen in figures 5 and 7, combining two 90° mitred ducts produces a 
rough approximation for an S-duct. As shown in figure 7, the rear and forward 
faces have been fitted with a steep S-curve. The modified double mitred duct 
has roughly the same characteristic as the single mitred duct, as seen in fig- 
ure 7. The energy reflection coefficient approaches unity near the first cut- 
off (but slightly less), at low frequencies the wave is passed unattenuated, 
and at higher frequencies the reflection coefficient is near unity. In addi- 
tion, a low frequency hump appears at a frequency of « near 1.5. Perhaps 
the duct is now acting analogously to an acoustic expansion chamber at low 
frequency which has similar types of humps due to resonance effects (ref. 12, 

p. 8). 


Example 3: S-Duct 

In the third example, the energy reflection coefficient was determined 
for the S-duct shown in figure 8 for various lengths. Generally, as the 
length increases the reflected energy decreases. In contrast to the mitred 
duct, the reflected energy approaches a minimum at the first node cut-on 
frequency of 3.14. Again, at the lower frequencies, the energy is nearly all 
transferred down the duct. The energy flux was determined by calculating the 
Poynting vector in the inlet and outlet duct from the calculated modal 
coefficients. 
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Example 4: Wall Absorbers 


Figure 9 shows duct configurations with an absorbing region along the 
upper and lower walls of the duct. Figure 9(a) show the straight wall config- 
uration while figure 9(b) shows the curved wall configuration. For the same 
absorbing material, the attenuation of both ducts will be compared. The wall 
material will have a relative permeability of 4.1 and a complex permittivity of 
1 - j 2 . 83 . An incident plane wave of dimensionless frequency co of 2ir is 
assumed. Figure 10 displays the finite-element grid used for both 
conf i gurations . 

For the straight duct, the axial component of the energy level (Poynting 
vector) is shown in figure 11. The drop in level is due to grazing absorption 
by the walls. Figure 12 displays contour plots of the magnitude of the mag- 
netic field for the straight duct absorber configuration. The solid line draw- 
ing shows the complete contours, while the symbols are plotted at selective 
points in the companion figure so that the actual magnitude of the contour can 
be readily determined. For convenience the magnitude of the magnetic field has 
been renormalized between 0 and 1 according to the formula 

IHI - |H | 


| x contour 

where H xmax is the maximum value of magnitude of the magnetic field in the 
plotting domain and H xm i n is the minimum value of H x in the plotting 
domain. 

As seen in figure 12, the magnetic field has decreased to near zero at 
the walls but remains relatively high in the center of the duct. In effect, 
electromagnetic energy beams through the center of the duct. This should not 
occur in the curved duct, because a more normal type impingement absorption 
wi 1 1 take place. 

For the curved duct, the axial component of the energy level is shown in 
figure 13. For the straight duct, the decrease in duct exit power due to wall 
absorption was -3.73 dB, while for the curved duct the decrease in exit energy 
due to wall absorption and reflection was -18.5 dB. The reflected energy at 
the inlet due to the change in duct curvature and wall impedance was -13.4 dB. 
The drop in energy level at the exit of the curved duct is considerably larger 
than the straight duct due to the more normal incidence of the magnetic field 
on the absorbing walls. Figure 14 displays the magnetic field contours. In 
this case, the curved wall prevents the beaming down the center which occurs in 
a straight duct. Consequently, the magnetic field quickly damps to near zero. 


xmax 


- H 


xmi n 


(32) 


CONCLUDING REMARKS 

A finite-element Galerkin formulation was developed to study transverse 
magnetic (TM) wave propagation in two-dimensional curved S-shaped ducts with 
both perfectly conducting and absorbing walls. The derivations from Maxwell's 
equations assumed that the material properties could vary with position result- 
ing in a nonhomogeneous variable property two-dimensional wave equation. This 
eliminated the necessity of finding the boundary conditions between the dif- 
ferent materials. Consequently, a complex structure can be easily modeled 
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simply by changing the property of elements in the finite-element domain. The 
reflection and transmission at the entrance and the exit of the curved duct are 
determined exactly by coupling the finite-element solutions in the curved duct 
to the eigenfunctions of an infinite, uniform, perfectly conducting duct. 

The numerical formulation is relatively simple to use and appears to give 
very accurate results. Example solutions are presented for a doubled mitred 
and S-duct of various lengths and with perfectly conducting and absorbing duct 
walls. The length of the S-duct is found to significantly affect the reflec- 
tive characteristics of the duct. Also, wall curvature is shown to greatly 
enhance the absorption properties of a duct. 
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AXIAL CORRDINATE, Z 


FIGURE 4. - A COMPARISON OF THE MAGNITUDE OF THE AXIAL 
MAGNETIC INTENSITY VARIATION ALONG THE LOWER WALL IN 
A UNIFORM DUCT WITH PERFECTLY CONDUCTING WALLS AND A 
CHANGE IN PROPERTIES OF THE MEDIA AT Z OF 0.25 (E f1 = 
1.0 AND e r2 = R.O) AS OBTAINED BY USING AN EXACT SOLU- 
TION AND A FINITE ELEMENT SOLUTION FOR A PLANE WAVE 
(MODE-ONE) INCIDENT AT Z = 0 WITH U = 27T. 



SINGLE 90° DOUBLE 90° MITRE S-DUCT 

MITRE 

FIGURE 5. - DOUBLE MITRE APPROXIMATION TO S-DUCT. 
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(B) CURVED DUCT. 

FIGURE 9. - ABSORBING WALL DUCT CONFIGURATIONS. 







(A) STRAIGHT DUCT. 



(B) CURVED DUCT. 

FIGURE 10. - DISCRETIZATION OF AIR FILLED WAVE GUIDE WITH 
ABSORBERS MOUNTED ALONG BOTH UPPER AND LOWER WALLS. 




-I.o -.5 0 .5 1.0 1.5 2.0 


AXIAL COORDINATE, z 

FIGURE 11. - EFFECT OF UPPER AND LOWER WALL ABSORBERS 
ON THE MAGNITUDE OF THE AXIAL FLUX OF ENERGY (POYNT- 
ING VECTOR) FOR A FIVE-MODE MODAL EXPANSION IN THE 
ENTRANCE AND EXIT DUCTS. 


RELATIVE 

MAGNITUDE 



FIGURE 12. - CONTOUR PLOTS OF MAGNETIC FIELD AMPLITUDE. 
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AXIAL COORDINATE, X 


FIGURE 13. - EFFECT OF UPPER AND LOWER WALL ABSORBERS ON 
THE MAGNITUDE OF THE AXIAL FLUX OF ENERGY (POYNTING 
VECTOR) FOR A FIVE-MODE MODAL EXPANSION IN THE ENTRANCE 
AND EXIT DUCTS. 
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